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Abstract. Dispersion effects induce new instabilities and dynamics in the weakly nonlinear 
description of light propagation in fiber Bragg gratings. A new family of dispersive localized pulses 
that propagate with the group velocity is numerically found and its stability is also analyzed. The 
unavoidable different asymptotic order of transport and dispersion effects plays a crucial role in the 
determination of these localized states. These results are also interesting from the point of view of 
general pattern formation since this asymptotic imbalance is a generic situation in any transport 
dominated (i.e., nonzero group velocity) spatially extended system. 

1. Introduction. Fiber Bragg gratings (FBG) are microstructured optical fibers 
that present a spatially periodic variation of the refractive index. The combination 
of the guiding properties of the periodic media with the Kerr nonlincarity of the fiber 
results in the very particular light propagation characteristic of these elements, which 
make them very promising for many technological applications that range from optical 
communications (wavelength division, dispersion management, optical buffers and 
storing devices, etc.) to fiber sensing (structural stress measure in aircraft components 
and buildings, temperature change detection, etc.), see, e.g., the recent review |5]. 

The amplitude equations that are commonly used in the literature to model one 
dimensional light propagation in a FBG are the so-called nonlinear coupled mode 
equations (NLCME) [17j [6j [7] , which, conveniently scaled,can be written as 

A+ - A+ = inA- + L4+(<r|^+| 2 + \A-\ 2 ), (1.1) 
A- +A~ = i K A+ +iA-(a\A-\ 2 + \A+\ 2 ), (1.2) 

where A are the envelopes of the two counterpropagating wavetrains that resonate 
with the grating, k is the strength of the coupling effect produced by the grating and 
a > is ratio of the self to cross nonlinear interaction coefficient (a = \ for a cubic 
Kerr nonlinearity [5]). The NLCME can be obtained from the full Maxwell-Lorentz 
equations using multiple scales techniques in the limit of small grating depth, small 
light intensity and slow spatial and temporal dependence of the field envelopes (see 
[7] for a detailed description of this derivation process) . 

It has been recently shown [TTJ [12] that light propagation in FBG can develop 
dispersive structures that are not accounted for in the NLCME formulation, and that, 
to correctly describe the weakly nonlinear dynamics of the system, the NLCME have 
to be completed with material dispersion terms: 

A+ - At = inA- + iA+(a\A+\ 2 + \A-\ 2 ) + \eA+ x , (1.3) 
A t - +A~ = inA+ + iA-(a\A-\ 2 + \A+\ 2 ) + \eA~ x . (1.4) 

The dispersive nonlinear coupled mode equations above (NLCMEd) are scaled as the 
NLCME: the characteristic length scale is the slow scale that results from the balance 
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of the advection term with the small effect of the grating, the characteristic time is 
the corresponding transport time scale (which sets to one the scaled group velocity), 
and the characteristic size of the wavetrains is the resulting one from the saturation of 
the small nonlinear terms. The slow envelope assumption forces the dispersive terms 
to be always small as compared with the advection terms; in other words, second 
derivatives of the slow amplitudes are much smaller than first derivatives. In the 
scaled equations this effect is contained in the scaled dispersion coefficient e (which 
measures the dispersion to transport ratio) and therefore, in order to be consistent 
with the slow envelope assumption, the NLCMEd must be considered in the limit 

The NLCMEd can be somehow regarded as asymptotically nonuniform, in the 
sense that the NLCMEd is an asymptotic model obtained in the e — ► (weakly 
nonlinear, slow envelope) limit that still contains the small parameter e. This is 
the unavoidable consequence of simultaneously considering two balances of different 
asymptotic order: one induced by the dominant effect of the transport at the group 
velocity (balance described by NLCME) and a second one that is associated with 
the underlying dispersive, nonlinear Schri'^i dinger- like dynamics of the system. The 
small dispersive terms in the NLCMEd are essential to describe the dynamics of the 
system when it develops small dispersive scales <5disp ~ y[ef. As it was shown in 
the NLCMEd in the e — > limit constitute a singular perturbation problem (cf. [10 ) 
and the onset of the dispersive scales is not a higher order, longer time effect; it takes 
place in the same timescale of the NLCME, no matter how small is the dispersion 
coefficient e. 

A solution of the NLCMEd for e = — 1CP 3 that exhibits small dispersive scales all 
over the domain is represented in Figure [TTT] note that, for short time, the dispersive 
structures just propagate with the group velocity but, for t ~ 1, they also interact 
with each other giving rise to a very complicated spatio-temporal pattern. The initial 
condition used in this simulation was a uniform modulus solution with a small random 
perturbation that, according to the dispersion-less NLCME formulation, was a stable 
solution. In order ensure that the small scales in Figure |1.1| are dispersion induced 
scales we have repeated the NLCMEd simulation but with a reduced dispersion s = 



10 /4. The result is plotted in Figure 1.2 where it can be seen that the small scales 



fill again the entire domain, but its typical size, <5di sp ~ y |e|, is now approximately 
one half of that in Figure |1.1| (see also the corresponding animations moviel . 1 .gif 
and moviel . 2 .gif). 

The main goal of this paper is to show that, in addition to complex spatio- 
temporal patterns, the dispersion effects can also give rise to new, purely dispersive 
localized states, which might be of interest from the optical communications point of 
view. It is interesting to note that the results in this paper apply also to Bose-Einstein 
condensates in optical lattices (a system that has recently received very much attention 
[151 118j ) and, in general, to any dissipation-less propagative system, extended in one 
spatial direction, reflection and translation invariant, and with a small superimposed 
spatial periodic modulation of its background, since the NLCMEd are the appropriate 
envelope equations for the description of the weakly nonlinear resonant dynamics of 
this kind of systems. 

In order to show that light propagation in a FBG can happen in the form of 
dispersive pulses, we derive and solve numerically in section 2 an asymptotic equation 
for a family of symmetric pulses, and, in section 3, we perform some numerical inte- 
grations of the complete NLCMEd to show that some of the pulses in this family do 
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Fig. 1.1. Space-time representation of a solution of the NLCMEd exhibiting small dispersive 
scales all over the domain (a = 1/2, k = 2, e = — 10 -3 , and periodicity boundary conditions). See 
the file moviel . 1 . gif for an animation of the onset of the dispersive scales. 
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Fig. 1.2. Space-time representation of a solution of the NLCMEd exhibiting small dispersive 
scales all over the domain (a = 1/2, ft = 2, e = — 10 — 3 /4, and periodicity boundary conditions). 
See the file moviel .2. gif for an animation of the onset of the dispersive scales. 
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Fig. 2.1. Sketch of a dispersive pulse on top of a continuous wave. 



propagate as stable localized structures. Finally, some concluding remarks are drawn 
in section 4. 

2. Dispersive pulses. The starting point is the continuous wave (CW) family 



of constant uniform modulus solutions of the NLCME ( 1.1 1-( 1.2 1 



A+ W = p cos 0e iut+imx , (2.1) 
Ac W = psin(? e - t+ ""* , (2.2) 



where p > is the light intensity flowing through the fiber, 9 G [—§,§] measures 
the relative amount of both wavetrains, and the frequency and wavenumber of the 
amplitudes are given by 

K a + 1 o , 

p 2 , (2.3) 



sin 26 2 

m = (^-- a -±p 2 )cos26. (2.4) 
sin 20 2 

The CW with \uo\ ~ 1 and |m| ~ 1 are approximate solutions of the NLCMEd (up 
to order e corrections) and its stability was first analyzed in [5] and then completed 
in [TT], where it was found that, for both signs of the dispersion coefficient, there 
are dispersively unstable CW which are stable in the dispersion less context of the 
NLCME. 

We now look for localized dispersive pulses propagating on top of one of the 



amplitudes of a stable CW, as sketched in Figure 2.1 In order to turn the background 
CW into a constant, it is convenient to first perform in the NLCMEd the following 
change of variables 



A+ 

' ~ e" 



to obtain 



- F+ + i(w - m)F+ = ikF~ + iF + (<j\F + \ 2 + \F~\ 2 ) +isF+ x , (2.5) 
F t ~ + F~ + i(w + m)F- - i K F+ + iF-{a\F-\ 2 + \F+\ 2 ) + ieF~ x . (2.6) 



A localized dispersive pulse on F + depends on the fast spatial scale X = x/yj\e\ and, 
according to (2.5)-(2.6), in the short time scale, T = t/\/\e\ ~ 1, it just propagates 



with the group velocity, suggesting that we have to look for solutions of the form 



F~ 



with r\ = X + T. Inserting the above ansatz into (2.5l-(2.6l yields 



F 0v = 0, that gives F = psin9, 

which means that F~ remains in first approximation equal to the unperturbed CW. 
Similarly, the following equation is obtained for Fq~ 



F+ + i(cu - m)F+ = mpsm9 + iF+(a\F+\ 2 + p 2 sin 2 9) ± iF ( 



07773 ' 



(2.7) 



where the + (— ) sign corresponds to e positive (negative), together with the boundary 
conditions 



Fq — ► pcos9 for 77 — > ±00, 



(2.8) 



which ensure that the background CW is recovered away from the pulse. 

For dispersive puls es propagating over a zero background we have to set to zero 
p, u) and to in eq. (2.7). A standard nonlinear Schr'i,;^ dinger (NLS) equation is then 
obtained, which is known to exhibit localized pulses (solitons) in the focusing case 
of positive dispersion (recall that a > 0). Note that in this case the effect of the 
grating is completely gone because it takes place only through the background state. 
The only nonzero Fourier spectrum components of these NLS solitons correspond, in 
first approximation, to very large dispersive wavenumbers (~ 3> 1) that are 

so off-resonanc e th at they simply do not feel the grating. 

Equation (2.7 1 is a nonlinear Schr'i^ ^dinger equation with a direct forcing term 



coming from the effect of the grating. The steady solutions of this equation and their 
stability properties were analyzed in (2J |3], where explicit analytic expressions for 
the steady pulses were found. In this paper we consider the more general family of 
traveling localized solutions. More precisely, we look for traveling pulses of the form 

F + = pcos6»(l + a(rj + vt)), 

where v represents, in the original variables, a small correction of the group velocity. 

br a, after making use of relation 
vt)(y/<rp cos 9), can be written as 



The resulting boundary value problem for a, after making use of relations (2.3 1 and 

(v- 



(2.4) and the rescaled variable £ 



a i6 



aa — ( I a. 1 2 + a + a) (I + a) , 

-> as £ — > ±00 



(2.9) 
(2.10) 



where v — v/(y / ap cos 9) and a = k tan(0) / {up 2 cos 2 9). We are restricting our search 
to the focusing case of positive sign in eq. ( |2.7| and, in order to have a dispersively 
stable background CW, we have to consider only the range < 9 < ? (see [TT]). 



which implies that we have to look for pulses in eqs. ( 2.9 )-( 2.10 1 only for a > 0. On 
the other hand, eqs. (2.9 1-( 2. 10 ) remain invariant under the transformation v — > — v 
and a 



a, and therefore we can also set v > 0. 
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Fig. 2.2. Eigenvalues of the trivial state as a function of the parameters a and v 2 . Sketches 
show the distribution of the four eigenvalues in each of the distinct regions. 



The solutions of eqs. (2.9 1-( 2. 10 1 correspond to traveling dispersive pulses and 



can be regarded as homoclinic orbits (in the variable £) connecting the trivial state 
(i.e., a — 0) back to itself. 

To analyze the existence of such connections we first consider the linearized system 
around the zero state: 



du 



(2.11) 



where u T = (a r , <2j, (<v)f > a r and a% are the real and imaginary part of a, 

{flr)i = da r /d^ and (a^)j = dai/d^, and the matrix Aac is given by 



1 \ 

1 

a-2 0-0 

a -v ) 



(2.12) 



The eigenvalues of this system are given by = ^fq± and Aj. = 

2a - (v 2 + 2) ± y/(v 2 + 2) 2 - 4v 2 a\ /2. 



rj±, with 



(2.13) 



Fig. 2.2 shows the behaviour of the four eigenvalues A]j_ and A^ in the (a, v 2 ) plane. 
There are four distinct regions, separated by the boundaries a = 2 and a — (v 2 + 
2) 2 /Av 2 : (a) four real eigenvalues, (b) two pairs of complex conjugate eigenvalues, (c) 
four purely imaginary eigenvalues and (d) two real eigenvalues along with two purely 
imaginary eigenvalues. In regions (a) and (b) the unstable and stable manifolds of 
the origin are two dimensional, while the equilibrium is nonhyperbolic in regions 
(c) and (d) where the center manifold is, respectively, four- and two-dimensional. 
Homoclinic orbits belong to both the stable and unstable manifold of the origin. We 
investigate below the presence of homoclinic solutions in cases (a) and (b), where these 
correspond to the intersections of two dimensional stable and unstable manifolds in a 
four dimensional space |16j . 



The problem (2.9 1-( 2. 10 > is invariant under the symmetry 



(2.14) 



that comes from the time reversing (Hamiltonian) and spatial reflection symmetries 
of the NLCMEd. We further restrict our search for dispersive pulses to the case of 
reflection-symmetric pulses, i.e., to pulses that satisfy: 



a(O = a(-0. 

If we now set the symmetry axis to £ = 0, we can reduce the problem to a semi-infinite 
interval £ € [0, +00 [ with the boundary conditions 

(a r ,ai) = (a ,0) and ((a r )f, (a 4 ) ? ) = (0, b ), at £ = 0, (2.15) 
(a r , di) — ► (Q, 0) as £^00. (2.16) 

Finally, to numerically compute the profiles of the symmetric pulses, we replace 
the infinite interval by a finite one, [0,L]. Following jS], the resulting boundary 
conditions at £ = L are obtained by requiring that the solution projects only onto 
the subspace spanned by the eigenvectors associated with the decaying eigenvalues 



of the matrix A^ (see eq. (2.11)). In summary, the boundary value problem that 



we integrate numerically is given by eq. (2.9 1, which we rewrite as a real first order 



system of four equations in ]0, L[ together with the four boundary conditions 

C* u = at £ = 0, and (2.17) 

CooU = at £ = L, (2.18) 

where u T = (a r , a i; (a r ) 5 , (a^), 

*-(!!!!!). <-) 

and Coo is a matrix whose rows are the left eigenvectors of associated with the 
exponentially growing directions 

c = ( (v+-a~v 2 ) <w/\A/+ (v+-a)/y/v+ v\ , 2 20) 
00 I (rj_ — a — v 2 ) av/yfrjZ (r]- — a) / \JfjZ v J ' 



For each value of a this problem is solved using a shooting method. We start from 
the known solutions for v = obtained in [3] and apply numerical continuation tech- 
niques to locate the propagating pulses with v > 0. This procedure for setting the 
boundary conditions at £ = L, rather than simply imposing a(L) — 0, allows the 
shooting method to converge faster, and the results obtained are found to be essen- 
tially independent of L for L > 10. 

The left panel of Fig. |2.3| shows several families of homoclinic orbits represented 
in the (a ,v) plane for different values of a and corresponding to the case a > 0. 
The dashed line separates the regions where the homoclinic orbit connects to a saddle 
point, and where it connects to a saddle- focus, while the open circles correspond to 



the solutions shown on the right panels. In Fig. 2.3 II a-c the solutions can be seen 



to develop oscillations as we move towards qq = 0. This corresponds to moving 



along a horizontal line in Fig. 2.2 (constant a) and to the right, approaching region 
(c) (precisely at ao = 0, v — ^/a + y/a — 2) where the eigenvalues become purely 
imaginary. 

This oscillatory behaviour, however, is not observed for the families of pulse so- 



lutions found for ap < 0, as seen in Fig. 2.4 Instead, these curves display turning 
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Fig. 2.3. (%) The solid lines correspond to homoclinic cycles to the origin with oo > for the 
indicated values of a. The dashed line bounds the regions where the origin is a saddle node (on the 
left) and saddle-focus (on the right). (II) Spatial profiles of three pulses for the values marked in (I) 
with open circles. 




points and, to better appreciate these limit points, the results have been plotted in 
the plane (E,v), where E is the positive quantity 



In this case, as we move along the curves for fixed a and past the turning point, the 



pulses develop two extra humps that tend to move away from the origin (see Fig. 2.4 
II (a) and (b)). 
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Fig. 3.1. Space-time representation of the solution of the NLCMEd (a = 1/2, re = 1, e = 10~ 5 , 
and periodicity boundary conditions) for a unstable pulse propagating on top of a CW (p = 1 and 
8 = tt/4). The pulse parameters are a = 4 and v = 1, and it corresponds to the point labeled (a) in 
Fig. \2.4\ See the file movie3. 1 . gif for an animation of this pulse destabilization. 



3. NLCMEd simulations. After having found a two parameter (a, v) family 
of symmetric dispersive pulses (DP) that can be numerically continuated from the 
solutions for v — obtained in [3], we now proceed to study their stability. 

The idea is not to perform a complete stability analysis of the family of DP 
but to show that stable DP can be found, and that the DP can thus be considered 
as robust realizable localized structures of light propagation in FBG. To do this we 
select several DP, place them on top of their corresponding background CW, add 
a small random perturbation, and use them as initial conditions for the full system 
of NLCMEd ( fOfLlI that we numerically integrate for a certain amount of time 
with periodic boundary conditions. The numerical method for the integration of the 
NLCMEd uses a Fourier series in space with Np modes and a fourth order Runge- 
Kutta scheme for the time integration of the resulting system of ordinary differential 

9 



equations for the Fourier coefficients. The stiff linear diagonal terms associated with 
the small dispersion coefficients are integrated implicitly, and the nonlinear terms are 
computed in physical space with the usual 2/3 rule to avoid the aliasing effects [I] 
(the maximum required resolution for the simulations in this paper were Np = 4096 
and At = .0005). 

Unstable DP simply do not persist and their shape changes, as it can be appreci- 
ated from the spatio-temporal evolution shown in Fig. |3 . 1 1 which corresponds to the 



pulse labeled (a) in Fig. 2.4 (see also the corresponding animation movie3 . 1 . gif ). 
Note that the size of the pulse grows, at t — 5.5 it is larger than at t = 3, and then 
it decays again at t — 7.0. A few time units later the oscillatory tails spread and 
the pulse structure is eventually lost (not shown in the Figure). All DP explored for 
a < propagated over the same simple C W with parameters p = 1 and 9 — ir/A (cf. 
eqs. (2.1 )-(2.2)), which corresponds to a = 4 in eq. (2.141, and all were found to be 



unstable regardless of its propagation speed v. 





Fig. 3.2. Space-time representation of the solution of the NLCMEd (a = 1/2, k = I, e = 10 — 5 , 
and periodicity boundary conditions) for a stable pulse propagating on top of a CW (p = 1 and 
9 = tt/4). The pulse parameters are a = 4 and v = 2.6, and it corresponds to the point labeled (a) 
in Fig. \2.3\ See the file movie3. 2. gif for an animation of this propagating pulse. 



On the other hand, for ao > and for the same background CW (i.e., a = 4), the 
DP are found to be unstable approximately for < v < 2.2 and stable for v > 2.2. 
The evolution of two stable pulses is shown in Figs. |3.2] and |3.3| which correspond to 
v — 2.6 and v — 3.2, respectively, where the structure of the slightly perturbed pulses 
is seen to remain now virtually unaltered after more than 40 time units (sec also the 
corresponding animations movie3 . 2 . gif and movie3.3.gif). 

Another very interesting feature of the DP that is worth mentioning is the fact 
that they are somehow transparent to each other: two DP propagating in opposite 
directions just pass through each other without distortion. This is illustrated in 



Fig. 3.4 (see also the animation movie3. 4.gif) where two stable DP (corresponding 
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Fig. 3.3. Space-time representation of the solution of the NLCMEd (a = 1/2, n = 1, e = 10 — 5 , 
and periodicity boundary conditions) for a stable pulse propagating on top of a CW (p = 1 and 
8 = The pulse parameters are a = 4 and v = 3.2, and it corresponds to the point labeled (b) 

in Fig. \2.3\ See the file movie3. 3. gif for an animation of this propagating pulse. 



to those in Figs. 3.2 and |3.3| ) are sent towards each other and after 40 time units 
(approximately 80 collisions) they still remain practically undistorted. The reason 
of this behavior is the dominant character of the transport effect induced by the 
group velocity If we rewrite the NLCMEd for a DP with short (dispersive) spatial 
and temporal scales X = x/y/e ~ 1 and T = t/y/e ~ 1, they take the form of two 
uncoupled wave equations 

A)p = A x + . . . , 

Arp — + . . . , 

and then it is clear that DP travelling in opposite directions simply propagate through 
different channels, and are (in first approximation) completely independent. 

4. Conclusions. In this paper we have studied the effect of dispersion in the 
weakly nonlinear dynamics of light propagation in a FBG. We have shown that the 
(often negelected) small dispersion terms play a crucial role in the transport dom- 
inated dynamics of light propagation in FBG. Dispersion can give rise to complex 
spatio-temporal chaotic states, but also to a new family of localized states (DP) that 
propagate over a CW background and do not interact with each other. These DP 
are approximately advected by the group velocity, but this transport effect does not 
play any role in the determination of their internal structure, which results basically 
from a balance between nonlincarity and dispersion. It is also important to emphasize 
that this type of dynamics is not contained in the standard dispersion-less NLCME 
formulation for light propagation in FBG. Moreover, this behavior is just the result of 
the competition of two effects with different asymptotic order: transport and disper- 
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Fig. 3.4. Space-time representation of the solution of the NLCMEd (a = 1/2, K = 1, s = 10~ 5 , 
and periodicity boundary conditions) showing the simultaneous propagation in opposite directions of 
the two stable pulses from Figs. \3.S\ and See the file movie3.4- gif for an animation of this 

pulse interaction. 

sion, and this is a generic situation that applies to any propagative extended system 
unless some special care is taken to reduce the group velocity (similar effects have 
been previously described in the context of Hopf bifurcation in dissipative systems 
1 3J and in parametrically forced surface waves [2]). 
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